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_0 I ABSTRACT Surface plasmon polaritons (SPPs) at the interface of a columnar thin film 

^ ' (CTF) and metal exist over a range of propagation directions relative to the morphology of 

(— I ■ the CTF which depends on the tilt of the columns in the CTF. The phase speed of the SPP 

CL|! wave varies mainly as a function of the tilt of the CTF columns. Both the confinement of the 
SPP wave to the interface and the decay of the SPP wave along the direction of propagation 

CN ' depend strongly on the direction of propagation relative to the morphologically significant plane 

, of the CTF. The greater the columnar tilt in relation to the interface, the shorter is the range of 

I propagation. Because of its porosity and the ability to engineer this biaxial dielectric material, 
the CTF-metal interface may be more attractive than traditional methods of producing SPPs. 
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1. INTRODUCTION 

The surface plasmon polariton (SPP) has been the object of intense investigation for several 
decades [1]. In recent years, the understanding of the phenomenon has been harnessed for a wide 
array of analytical and biomedical applications [21 [3l H] • Although various configurations are 
used to launch and detect the SPP, in the Kretschmann configuration, a metal film is interposed 
between a low-refractive-index dielectric medium and a high-refractive-index dielectric medium 
[5]. A light beam is launched in the high-refractive-index medium to impinge on the planar 
interface between that medium and the metal. The SPP travels along the planar interface of 
the metal and the low-refractive-index medium. 

The propagation of SPPs at the planar interface of a metal and an anisotropic dielectric 
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Figure 1: Geometry of CTF-metal interface and SPP-wave propagation. 

material has been examined theoretically [6]. For this communication, we investigated the 
propagation of the SPP along the interface of a metal and a columnar thin film (CTF), in relation 
to the morphology of the CTF. The columnar thin film, an artificial material, is effectively an 
optically biaxial continuum [7], and porous [SI Ch. 6]. The porous CTF may offer a medium in 
which to embed analyte and/or recognition molecules to which the analyte may bind. 

A columnar thin film is produced by directing a vapor at an angle Xv to a substrate. Under 
suitable conditions, parallel columns form at an angle X ^ Xd to the substrate; see Figure [TJ 
The columns are composed of multimolecular clusters with ~ 3 nm diameter which, in turn, 
form columns with ~ 100 nm cross-sectional diameter, depending on the evaporant species and 
the deposition conditions. The geometry considered here is the same as that in a predecessor 
paper [9], except that the substrate is now made of a metal instead of a transparent dielectric, 
and is also illustrated in Figure [TJ 

The plan of the communication is as follows. Section 2 contains a description of the geometry 
of the problem and the constitutive relations. For the details of the field representation and the 
derivation of the dispersion relation for the SPP wave, we refer the reader to the predecessor pa- 
per [9]. Numerical results for a particular CTF-metal interface, titanium oxide CTF-aluminum, 
are presented in Section 3. Concluding remarks are presented in Section 4. 

The following assumption and notational conventions are used. An exp(— iwt) time-dependence 
is implicit, with uj denoting the angular frequency. The free-space wavenumber, the free-space 
wavelength, and the intrinsic impedance of free space are denoted by ko = Wy^e^/Xo = ^/co-, 
Xo = It: jko^ and r/o = ^ ^o/^o-, respectively, with /io and eo being the permeability and permit- 
tivity of free space. Vectors are underlined, and dyadics are underlined twice. Cartesian unit 
vectors are identified and M,. 



2. GEOMETRY AND CONSTITUTIVE RELATIONS 

The geometry of the problem is illustrated in Figure [TJ With the interface between the metal 
substrate and the CTF denoted by z = 0, the metal occupies the half-space z < 0, while the 
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CTF occupies the half-space z > 0. Without loss of generality, the direction of propagation is 
taken to be parallel to the x-axis, and the angle between the direction of propagation and the 
morphologically significant plane [HJ Ch. 7] of the the CTF is denoted by ^ . 

The relative permittivity dyadic of the CTF may be stated as |8i Ch. 7] 

where na^b,c are the principal refractive indexes and the unit vectors 

Un — ~ (Mx cos tp + Uy sin tp) sin x + ILz cos X 1 

= (Uj, cos + My sin -0) cos X + sin X > . (2) 
= —Uy cos ip + Ux sin J 

All four quantities — nafi^c and the column inclination angle x — depend on the evaporant 
species and the vapor incidence angle Xv- The refractive index of the metal is denoted by n^. 

Expressions for the field equations and dispersion relation have been worked out elsewhere 
[9]. It suffices to mention that the wave vector in the substrate may be written as 

= ko - iqs Uz) ^ (3) 

where x and 

Qs = +\/x2 - (4) 

are the normalized propagation constant and the decay constant, respectively. In order for the 
SPP wave to be confined to the interface z = 0, we must have Re [qs] > 0. Unlike for the 
CTF-dielectric interface (9], x is expected to be complex valued since the dissipative properties 
of the metal can not be ignored. Similarly, in the half-space z > occupied by the CTF, the 
wavevector may be written as 

k^ = ko {xUx + '''QclLz) ) (5) 

and again we must have Re [qc] > for localization of the SPP wave to the interface. The SPP 
wave comprises two partial waves in the CTF, identified by qd and qc2- 



3. NUMERICAL RESULTS 

Hodgkinson et al. [10| measured the constitutive relations for CTFs made of oxides of tantalum, 
titanium, and zirconium. We chose titanium oxide to illustrate SPP propagation at a planar 
CTF-metal interface. The empirical relationships determined at Ao = 633 nm for the titanium- 
oxide CTF can be written as 

na = 1.0443 + 2.7394 [yj^ " 1-3697 {^^J^ > (6) 

rib = 1-6765 + 1.5649 - 0-7825 > (7) 

He = 1.3586 + 2.1109 ( - 1.0554 ( , (8) 
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Figure 2: Phase speed of the SPP wave at the planar interface of a titanium-oxide CTF and an 
aluminum substrate, relative to the speed of light in vacuum, versus ■0 for Xv = 7.2° , 20°, 60°, 
and 90°. 



and 

tanx = 2.8818 tan x^, (9) 

where Xv and x are in radian. We must caution that the foregoing expressions are applicable 
to CTFs produced by one particular experimental apparatus, but may have to be modified for 
CTFs produced by other researchers on different apparatuses; hence, we used these expressions 
for the numerical results presented in this section merely for illustration. Metals commonly used 
for SPP systems are aluminum, copper, silver, and gold. For illustration, we chose aluminum 
with a complex refractive index Ug = 1.38 + 7.61z at Aq = 633 nm [llj. We note that values of 
X for SPP propagation when = 0° have been reported elsewhere |12j . 



Calculations were carried out for Xv = 7.2°, 20°, 60°, and 90°. The lowest value Xv = 7.2° 
corresponds to x = 20° and represents the approximate, currently achievable lower limit of x- 
As previously stated, the wave number k = ko>c along the direction of propagation must be 
complex valued. The real part of k determines the phase speed v = a;/Re(K) of the SPP wave. 
The phase speed relative to the speed of light in vacuum is calculated as 

V 1 

Vrel = — = p . s • (10 

Co Re(xj 

Figure [2] shows function of for Xv = 7.2°, 20°, 60°, and 90°, with ip restricted to the 

range [0°,90°]. The curves are drawn over the -i/'-ranges for which the boundary conditions at 
z = could be satisfied and thus represent the ranges over which a SPP wave can exist. For 
each value of ip shown, at which a SPP wave exists, there also exist three other values at —■0, 
180° + 0, and 180° — ip with an identical value of Vrei- The relative phase speed Vrei decreases 
as Xv increases with the most rapid change occurring at the low end of the Xi;~range. 

The greatest variation of function of ip is observed for Xv = 20°. This curve has a 

downward slope as Vrei decreases from roughly 0.540 to 0.516 as ^/J changes from 0° to 52°. At 
Xv = 60°, a much smaller change is observed: Vrei decreases from roughly 0.4050 to 0.4045 as 
ip changes from 0° to 90°. In contrast, the curve for Xv = 7.2° shows an upward slope as Vrei 
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Figure 3: e-folding distance of the SPP wave along the x-direction, the direction of propagation, 
versus ^ for the same conditions as in Figure [2j 



changes from roughly 0.735 to 0.738 as ip varies from 0° to 15°. Of course, the curve for Xv = 90° 
is completely flat as required by symmetry since Eqs. (l6|)-([9]) predict that the generally biaxial 
CTF becomes uniaxial at Xv = 90° with the axis of symmetry perpendicular to the interface 
z = 0. 

The imaginary part of k determines the rate at which the SPP wave decays along the di- 
rection of propagation. The e-folding distance for decay of the wave amplitude relative to the 
wavelength Aq may be calculated as 

_ 1 

27rlm(x) 

Figure [3] shows function of ip over the restricted range [0°,90°]. As with Vreh for each 

value of ip displayed, an identical value of Xg/ is obtained at —ip, 180° + ip, and 180° — ^p. As 
Xv increases, x^f decreases and the curve ofxef vs. tp flattens. At Xv = 60°, the curve is nearly 
flat with only ~ 2% change over the entire ■0-range. In contrast, at Xv = 7.2° and 20°, X(,f 
changes by nearly two orders of magnitude. 

Figures [2] and [3] let us conclude that an increase in the vapor incidence angle (i) reduces 
the phase speed and (ii) increases the attenuation of the SPP wave. Thus, a high value of 
Xv is inimical to long-range propagation for all ip, an understanding previous obtained only 
for the case of the direction of propagation lying in the morphologically significant plane (i.e., 
■0 = 0°) [12]. The properties of the SPP wave along the z-axis, perpendicular to the direction 
of propagation, are described by 

(i) Qci and (ii) qc2 for the two partial waves in the CTF, and 
(iii) Qm in the metal. 

Similar to Xej, three e-folding distances in the z-direction relative to the wavelength in vacuum 
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Figure 4: e-folding distances of the SPP wave relative to Aq in the z-direction for the same 
conditions as in Figure [2j a) partial wave 1 in CTF; b) partial wave 2 in CTF; c) in aluminum. 



may be calculated as 

^e/l = 27rlm{gei) | 
^e/m = 27rlm(g„) J 

Figure S] shows these three quantities as functions of ip over the restricted range [0°, 90°]. These 
three quantities describe the localization of the SPP wave to the interface. In Figure U^, the 
curves of^eji are similar to those of Xe/ displayed in Figure [3l with the curves for both Xv = 60° 
and 90° essentially flat. At Xv = 7.2° and 20°, Zg/i varies by about two orders of magnitude. 
Thus, for the two values of Xv for which the ■0-range for surface wave propagation ends and 
does not continue to ^lJ > 90°, the first partial-wave component of the SPP wave in the CTF 
becomes delocalized as ip approaches the end of its range. 

In Figure Hb, the curves of versus tp are also nearly flat for both Xv = 60° and 90°. At 
Xv = 7.2° and 20°, Ze/2 decreases as ip increases. The second partial-wave component, then, 
becomes more localized as tp approaches the end of its range. 

The e-folding distance in the metal, ^e/mi is essentially flat for all values of Xv in Figure Hh. 
The greatest variation of Zg/m. occurs when Xv = 20°, but is still less than 0.4%. 

Unlike Vrei and x^j, the values of the e-folding distances in the z-direction for the CTF are 
not identical at ip, —ip, 180° + ip, and 180° — V'- The value of any e-folding distance in the z- 
direction at ip is, however, equal to the value at —ip, and the value at 180° +ip is equal to that at 
180° — tp. Since the difference in most cases is slight, we present the percent difference between 
the values at ip and —ip, and the values at 180° + ip and 180° — ip in Figure [5j The % differences 
for the e-folding distance of the first partial-wave component in the CTF, Zgfi, are shown in 
Figure The curve for Xv = 7.2° shows the largest difference of about 13% and occurs at 
the end of the t/^-range where ip = 15°. The difference decreases rapidly with increasing Xv At 
Xv = 60°, the maximum difference is only about 0.1%. Differences in e-folding distances for 
the partial-wave component 2 in the CTF are shown in Figure [5)3 with the maximum observed 
value of 1.5% — almost an order of magnitude smaller than for the partial-wave component 1. 
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Figure 5: The % difference between Zefi{ip) and 2;e/2(180° +ip) for same conditions as in Figure 



Differences in e-folding distances appear to decrease as Xv increases with a maximum difference 
Xv = 60° of only about 0.01%. However, the decrease is not monotonic as the maximum 
difference at Xv = 20° is about 50% larger than that at Xv = 7.2°. 



Examining SPP-wave propagation at a planar CTF-metal interface, with demonstrated empir- 
ical data to characterize both the CTF and the metal, we found that the direction, denoted by 
-0, at which SPP-wave propagation is possible is limited, in general, and depends on the tilt 
angle x of columns in the CTF. This tilt angle may be predicted from the vapor incidence 
angle Xv set during the manufacture of the film. The ■0-range for SPP waves increases as Xv 
increases. At sufficiently large values of Xv, SPP propagation becomes possible in all directions. 
Even at the lowest value {xv = 7.2°) examined by us, the ip~range is greater than 10°. This is 
in contrast to surface-wave propagation at the interface of a CTF and an isotropic dielectric 
material, for which the ■0~range is only a fraction of a degree for all values of x [9j . 

The phase speed of the SPP wave shows a strong dependence on Xv, but is also mildly 
dependent on ^p. The e-folding distance along the direction of propagation is strongly dependent 
on ■0 at low values of Xv , varying by several orders of magnitude, but is relatively fiat at larger 
values of Xv A high value of Xv is inimical to long-range propagation for all ip. At the end of 
the '(/'■range, for those values of Xv for which the propagation directions are limited, the SPP 
wave becomes delocalized from the interface on the CTF side. This is seen as an apparent 
divergence in the e-folding distance in the direction perpendicular to the interface for one of 
the two partial waves in the CTF as ip approaches the end of its range. The e-folding distance 
perpendicular to the interface for the other partial wave in the CTF, on the other hand, tends to 
decrease as ip approaches its limiting value. In the metal, the e-folding distance perpendicular 
to the interface does not vary much with ip. 

The demonstration of a SPP wave at the interface of a CTF and a metal offers certain options 
in SPP technology. With a wide choice of possible evaporant materials [13] and vapor incidence 
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4. CONCLUSION 
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angles, the CTF-metal interface offers a large latitude in the engineering of SPP systems. The 
porosity of the CTFs may also offer some advantages. With its nano~scale porous structure, 
a CTF could allow molecular-scale analytes to reach the interface for analysis while screening 
out larger particulates. The detection of viruses has become important in recent times. CTFs 
with pores on the nano-scale may allow viruses access to the interface while excluding larger 
organisms. Through photolithographic techniques it is possible to pattern CTFs [14J. Various 
samples could be placed, possibly in the field, on isolated CTF patches on a single substrate for 
later analysis. Also, different patches could be embedded with different recognition molecules 
which could be used to perform multiple searches on a single sample. Other options have yet 
to be imagined. 
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